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Spin and charge-current dynamics after ultrafast spin-polarized excitation in a normal metal are 
studied theoretically using a wave-diffusion theory, ft is shown analytically how this macroscopic 
approach correctly describes the ballistic and diffusive properties of spin and charge transport, but 
also applies to the intermediate regime between these two limits. Using the wave-diffusion equations 
we numerically analyze spin and charge dynamics after ultrafast excitation of spin polarized carriers 
in thin gold films. Assuming slightly spin-dependent momentum relaxation times, we find that a 
unified treatment of diffusive and ballistic transport yields robust signatures in the spin and charge 
dynamics, which are in qualitative agreement with recent experimental results [Phys. Rev. Lett. 107, 
076601 (2011)]. The influence of boundary effects on the temporal signatures of spin transport is 
also studied. 

PACS numbers: 72.25.Ba, 73.23.Ad, 75.78.Jp, 85.75.-d 



I. INTRODUCTION 

There has been a persistent interest in spin transport 
phenomena because of their importance in physics and 
its application in spintronic devices. Spin transport in 
magnetic multilayers composed of ferromagnetic metals 
and simple metalsOhas been perhaps the most intensely 
studied of these phenomenaP^ for instance, spin injec- 
tion/detection using currents flowing perpendicular to 
the plane of the layers (CPP)J^spin manipulation by elec- 
tric and magnetic fields^ and CPP giant magnetoresis- 
tance (GMR)PFor a macroscopic description of spin cur- 
rents and the accompanying spin accumulation at layer 
boundaries, one often applies macroscopic spin-diffusion 
theorypHHD which yields an unphysical infinite propaga- 
tion velocity for a spin signal. This property of the spin 
diffusion equation becomes troublesome at high frequen- 
cies and/or switching speeds, and can be circumvented 
by using a macroscopic spin wave- diffusion equation, as 
discussed in Refs.fTTIandlrH In this paper, we are mainly 
concerned with (i) a general discussion how the wave 
character of spin transport is related to ballistic trans- 
port in the macroscopic spin wave-diffusion theory and 
(ii) an application of the theory to the case of metallic 
heterostructures in which the spin polarization is created 
by ultrashort optical pulses. The motivation for the lat- 
ter problem comes mainly from recent experimental re- 
sults,^ where spin and charge resolved transport in thin 
gold films over 100 nm after ultrashort-pulse excitation 
was observed and signatures of both ballistic and diffu- 
sive transport were found. While the results on charge 
dynamics in Ref. UJ2 bear some resemblance with earlier 
results,^ the measured spin dynamics yield additional 
characteristic signatures, which we analyze in the frame- 
work of spin wave-diffusion equations. Since our theoret- 
ical approach contains ballistic and diffusive dynamics as 
special cases, it is well suited for a study of the behavior 
at the transition between ballistic and diffusive transport. 



The paper is organized as follows. The time-dependent 
equation system, which forms the basis of our discussion 
is reviewed in Sec. [TTJ In Sec. |III| we show analytically 
how the wave character of these equations is related to 
ballistic transport. In Sec. IV we present numerical re- 
sults for the spin and charge dynamics in a gold film 
after spin-polarized optical excitation, and compare our 
results to recent experiments. We summarize our results 
in Sec. El 



II. TIME-DEPENDENT EQUATION SYSTEM 

We base our discussion of the spin dynamics in a metal 
on the dynamical equations for the macroscopic spin- 
current density J s (z,t) and spin density n s (z,t), as de- 
rived in Ref. [TTJ 

dJ s (z,t) dn s (z,t) n s (z,t) - n- s (z,t) 

dz + dt r sf ' ( ' 

J s (z,t) _ _ 2 dn s (z,t) _ dJ s (z,t) 

t s ~ Cs dz dt • ( ] 

The equations describe macroscopic spin and charge cur- 
rents carried by conduction electrons in parabolic bands 
with spin s = ±|. For the description of multilayers it 
is assumed that the same bands with spin label s ex- 
ist in each part of the heterostructure. Further, r s f and 
t s are the spin- flip and the (spin-dependent) momentum 
relaxation times, respectively. Eq. ([I]) for the charge cur- 
rent J and charge density n has a similar form; only the 
flip term on the RHS of Eq. ([!]) is absent. The equation 
corresponding to Eq. |lj) involving the charge current and 
density is just the continuity equation. 

In the derivation of Eq. ([I]) from the Boltzmann trans- 
port equation, conductivities and momentum relaxation 
times are introduced to describe dynamical changes of 
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the electronic distributions close to the Fermi energy. In 
ferromagnets, the momentum relaxation times t s are in 
general different for majority and minority electrons, but 
in non-magnetic metals they are usually assumed to be 
spin-independent. To rigorously include excitation fur- 
ther away from the Fermi energy in Eqs. ([!]) and (|2|), one 
would have to augment them by dynamical equations for 
the electronic energy density (or temperature). We do 
not do this here, but observe that electrons excited at 
high energies will relax quickly towards the Fermi en- 
ergy, so that we can phenomenologically include effects 
of a spin- dependent excitation at energies further away 
from the Fermi energy by allowing for spin-dependent 
relaxation times in the non-magnetic metal. Regardless 
of the excitation conditions, a spin-dependence of the 
momentum relaxation-times may also arise in an experi- 
ment, such as the one in Ref. 1131 if an external field for 
the measurement of magneto-optical signals is applied. 

To make this paper more self-contained, we collect 
some results from Ref. [TTJ There it was shown that an 
important consequence of the macroscopic wave-diffusion 
equations is the occurrence of the spin-signal-propagation 
velocity, 
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where i>f is the Fermi velocity. Eqs. and |2]) can 
be combined to yield a wave-diffusion equation 11 for the 
total spin density n m = n + — rt_, 
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where Ty = ^ is the spin relaxation time and r the 
"average" momentum relaxation time. The equation 
above contains a second-order time derivative, which is 
absent in the spin diffusion equationP This additional 
term shows that it takes a finite time for the spin cur- 
rent toadjust to the gradient of the spin accumula- 
tion! 15 * 1 ^ The second-order time and space derivatives 
lead to a wave character of the dynamical spin and charge 
transport in addition to its diffusive character and that 
time-dependent spin transport should be understood at 
a macroscopic level in a wave-diffusion duality picture. 
From this viewpoint, the time-dependent spin diffusion 
equation can be regarded as an approximation of the 
wave-diffusion behavior in the long-time/low frequency 
limit. 



III. BALLISTIC AND DIFFUSIVE TRANSPORT 
FROM THE WAVE-DIFFUSION EQUATION 

The wave-diffusion equations described above yield 
time and space dependent distributions functions, 
n m (z, t), that will be discussed in detail in the next sec- 
tion. Here we wish to connect the wave-diffusion picture 
of spin transport with its diffusive or ballistic character. 



To this end, we analyze the spatio-temporal behavior of 
the spin polarized carrier distribution by calculating the 
mean-square displacement!^ 



A*(t) = 



1 



dzz 2 n m (z, t). 



(5) 



with N m — dz n m (z, 0). We will interpret A 2 {t) — 
A 2 Z (0) oc t 2 as ballistic and A 2 z {t) - A 2 Z (0) cx t as diffusive 
transport in accordance with electronic wave packet dy- 
namics. To simplify the discussion further, we specialize 
to the case of a spin density that is initially concentrated 
at z = 0, i.e., the initial condition 



n m (z,0) = N m S(z), 



(6) 



so that A 2 Z (0) = 0. 

Following Weissp^we construct a solution of the wave- 
diffusion equation with the initial condition ^ using 
damped dispersive solutions of the form 

n m oc *). (7) 

In Ref. IIR k(ui) = k r (u>) + iki(uj) was chosen as a com- 
plex quantity to obtain the damping length as function of 
frequency. Here we let uj{k) become a complex function 
of the wave vector to determine the temporal damping 
characteristics. Substituting Eq. |7]) into Eq. Q, we get 
the dispersion relation 
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which describe an exponential decay of n m (z,t) and 
J m (z,t) in time. The initial condition ^ can now be 
satisfied by the Fourier integral 



n m (z,t) = N„ 
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from which one reads off the Fourier transformed density 



h m (k,t) = N m dze lkz n m (z,t) 

J — OC 
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The mean-square displacement can be expressed with the 
help of the second derivative of the Fourier transforma- 
tion and explicitly evaluated as follows 

1 d 2 n m (k,t) 
N m 
2c 2 t 



dk 2 



k=0 
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FIG. 1. Mean-square displacement A z (t) of the spin density 
vs. time: Transition from ballistic to diffusive spin transport. 
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FIG. 2. Same as Fig. [l] on a longer time scale. 



Eq. ( 12 ) describes the behav- 



where we assumed t <T\. 
ior of A|(i) for all times. 

In Fig. [T] we plot the time-dependence of the mean- 
square di splace ment for the parameters (t, T\, c s ) of gold 
(see Sec. |IV A[ ) . On time scales for which the spin relax- 
ation is not important, i.e., t < T\, Eq. (12) simplifies 
to 



A^)~2c> 
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For short times, t <C r, we have A|(i) oc (? s t 2 . The char- 
acteristic ballistic behavior of A 2 is consistent with the 
limit t -C t, which describes a time scale shorter than 
typical (momentum) relaxation processes. For times on 
the order of the momentum relaxation time r, momen- 
tum scattering processes come into play, and the ballistic 
behavior is gradually lost. More precisely, for r < t < T\, 
the behavior of the second moment is A 2 z (t) oc c 2 rt 
with corrections of the order exp(— t/r). By virtue of 
c 2 t = D and disregarding exponentially small correc- 
tions, this means A^(t) oc Dt, which is a signature of dif- 
fusive transport with diffusion constant D. Finally, for 
times on the order of (and longer than) the spin relax- 
ation time T\ we have an exponential decay of n rn (z,t) 
and also an exponential decay of A 2 (t). This result is 
simply a consequence of the relaxation of the spin polar- 




FIG. 3. (color online) Schematic setup of the optical pump- 
probe experiment on a Fe/Au heterostructure. 



ization due to spin-flip processes. These considerations 
clarify how the macroscopic wave-diffusion equations 
and ([2]) are capable of describing both the ballistic and 
diffusive properties of spin transport, as well as the tran- 
sition between the two limits. We will use these concepts 
in the discussion of the numerical results for spin trans- 
port after ultrashort optical excitation in the following 
section. 



IV. NUMERICAL RESULTS 

A. Model and computational setup 

In this section, the time-dependent equation system 
presented in Sec.[ll]is applied to a model setup that aims 
at analyzing the spin and charge transport as it occurs in 
the experiment of Melnikov et al.f^ as shown in Fig. [3] 
In the experiment, spin polarized carriers are excited op- 
tically by a 35 fs pump pulse in a Fe/Au slab structure. 
The excited carriers propagate from the iron into the ad- 
jacent nonmagnetic gold slab, and the charge and spin 
dynamics at the right gold surface are extracted from 
the time- dependent magneto-optical response at the gold 
surface. 

We are interested in the space- and time-resolved spin 
and charge dynamics in the gold layer, and introduce 
the following simplifications. First, we assume that the 
spin signal excited in the iron layer can be modeled as 
a boundary value for the spin and charge densities in 
the adjacent gold slab. Second, we neglect the spin 
and charge transport of holes. We then use the macro- 
scopic equations for the electron densities, n s (z,t), and 
current densities, J s {z,t), in the spin-up and spin-down 
bands, respectively, of a ferromagnet-metal heterostruc- 
ture. These dynamical macroscopic equations can be de- 
rived from the Boltzmann transport equation for the mi- 
croscopic carrier distributions under the assumption that 
the deviations from an equilibrium distribution around 
the Fermi surface are smallP^ We therefore use the fol- 
lowing simplified view of the complicated injection and 
relaxation process: Spin-polarized carriers are created in 
the iron layer and move into the spin-degenerate gold 
bands where this is allowed in the band structure. We as- 
sume a very fast initial electronic energy equilibration, so 
that we have quasi-equilibrium distributions in the spin 
degenerate bands in gold, but with different densities. 



4 



The momentum relaxation-times of the quasi-equilibrium 
electrons of different spin, r + and r_, which arise micro- 
scopically from averages over the energy dependence of 
the electronic distributions, should then also be different. 
We therefore treat the difference of the relaxation times, 
r+ and r_, as a parameter that needs to be specified. In 
addition, we could also treat the spin polarization of the 
injected electrons as a parameter, but we choose not to 
do so to avoid the introduction of fit parameters. Since it 
turns out that the dependence of our results on the dif- 
ference of momentum relaxation-times is extremely weak, 
the only important parameters of our macroscopic equa- 
tions have values that are well established from electri- 
cal transport measurements: the Fermi velocity wp> the 
momentum relaxation time r s and the spin relaxation 
time T\. Using our approach thus yields the dynamics 
due to quasi-equilibrium transport, and can be compared 
to more involved treatments of the relaxation processes 
away from the Fermi energy. For instance, these relax- 
ation processes may lead to super-diffusive transport.^ 
In order to solve Eqs. and |2]) numerically for elec- 
tronic transport in the gold layer, we need the boundary 
condition at the interface to the ferromagnet, i.e. the left 
boundary of our computational domain, which we take to 
be z = 0. For simplicity, we assume that the spin-current 
density as a function of time is known at the Fe/Au in- 
terface. This approximation neglects the spin current 
flowing back to the iron layer due to the spin accumula- 
tion in the gold layer. At the Fe/Au interface, z = 0, the 
spin-current density is assumed to be a Gaussian 

J s (z = 0,t) = J°e- t2 /*«, (14) 

where J° is the maximum current density. Note that 
we assume the same initial condition for the currents in 
the minority and majority channel. Further, to is re- 
lated to the full-width at half maximum of the pulse via 
to = iFWHM/(21og2). We choose £fwhm = 35 fs and 
J° = 10~ 4 nm~ 2 fs . For the Au layer, we choose a mo- 
mentum relaxation time of r p = 30 fs (from Ref. 120"]) 
and c s = 0.808 nm/fs (v F = 1.4 nm/fs, see Re f. [Mil. 
The spin relaxation time Ti = 515 fs is determine d 6 * 11 * 21 * 
from the spin-diffusion length Z s f = 100 nm via the iden- 
tity Z s f = c s \J tT\ . The resulting electron mean- free path 
is A = v-pTp = 42 nm. 

In the following, we take the literature value of the 
momentum relaxation-time for one spin channel, i.e., 
t + = t p , but t_ 7^ t + to mimic the effect of different mi- 
croscopic distribution for majority and minority electrons 
that propagate into the gold layer. Some justification for 
our approximation is provided by the weak dependence 
of the numerical results on the actual value chosen for 

T_. 

We study two simplified variations of the setup de- 
picted in Fig. [3j In the first case, we assume an "infinitely 
thick" gold layer to avoid the influence of the reflected 
spin signal. In this case we assume that the thickness L 
of the gold slab is so large that spin and charge signals 
have not propagated to the right boundary during the 
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FIG. 4. Time evolution of the density n+(z,t) of majority 
electrons (a) and the total spin density n m (z,t) (b) at posi- 
tions z = 25 nm (solid), 50 nm (dashed), 100 nm (dotted), and 
150 nm (dash-dotted) of an infinitely thick layer. The excita- 
tion current pulse J s is also shown in arbitrary units(thin red 
line). 



time interval of interest, i.e., L l s {. In this time frame, 
the choice of the boundary condition at the right of the 
gold layer does not influence the results. For mathemat- 
ical convenience, we require both spin-(charge-) current 
density and spin (charge) density be zero. We also study 
the case of a finite thickness L of the gold layer. Here, 
we require the spin and charge current densities to van- 
ish at the right boundary of the gold layer. Although 
it may seem that this boundary condition applies to the 
experiment best, it leads to multiple reflections that do 
not show up in the experiment, as is discussed in more 
detail below. 

The numerical solution of Eqs. and ([2| employs the 
method of characteristics as detailed in Ref. ITT1 



B. Spin and charge dynamics 

We first discuss results for the time evolution of the 
spin and charge density in an "infinitely thick" layer. 
Fig. |4](a) shows the spin density n + (z = L,t) for major- 
ity electrons at different positions z ~ L, ranging from 
L = 25 nm to L = 150 nm. As mentioned above, the 
results at position z — L mimic the dynamics at the 
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surface for a slab thickness L. We find a delayed sharp 
rise and a peak in the carrier density. We hrst discuss 
the curve for z = 25 nm in Fig. |4]Ja), which is in the 
regime where ballistic transport properties are dominant, 
as shown by Fig. [I] The pronounced peak occurs at time 
ip Ca k = L/c = 34 fs, and therefore is related to wave- 
like propagation dynamics, which corresponds to ballis- 
tic transport as shown in Sec. |III| Further, the wave-like 
dynamics of the carrier distribution preserves to some ex- 
tent of the Gaussian shape of the driving current pulse 
J s (z = 0, t), which is created by the drive current at 
the left boundary. The difference between the full and 
the purely ballistic result is, of course, due to momen- 
tum and spin-flip scattering. These contributions lead 
to an asymmetric peak and the slow decay, but not to a 
broadening of the peak. This can be seen by comparing 
with the result for an infinitely thin layer (not shown) , in 
which the dynamics are completely ballistic. For peaks 
at larger distances (and longer times) the peak-like sig- 
nature is gradually lost as one reaches the regime where 
the diffusive behavior becomes dominant and the ballis- 
tic (wave-like) component is essentially only seen in the 
delayed onset of the signal. 

We therefore interpret the dynamics in Fig. [4](a) as 
ballistic for short times, as more or less diffusive for "in- 
termediate" times, and eventually as an exponential de- 
cay for long times. From the ballistic motion of electrons 
traveling with v-p in a straight line towards the opposite 
surface one may naively expect the arrival of a signal 
with the Fermi velocity instead of c s — Vp/y3. This 
is only true for carriers "going straight" to the opposite 
surface, and these carriers are responsible for the onset of 
the signal. The occurrence of the peak with the smaller 
velocity of c as well as the finite width of the peak shows 
that there is a spread in the time the carriers need to 
traverse a fixed distance in z direction. 

The curve for minority electrons n_ (z, t) (not shown) is 
almost indistinguishable from the one shown for n+(z, t) 
on the scale of Fig. |4|a) if we choose the same exci- 
tation, i.e., boundary condition, and change the mo- 
mentum relaxation time for "— " electrons by 5% to 
t_ = 31.5 fs. However, as Fig. |4^b) shows, the total 
spin density n m (z,t) = n + (z,t) — n_(z,t) is influenced 
quite drastically by this change. Qualitatively, this can 
be understood by considering the following picture of the 
dynamics of individual electrons, even though the calcu- 
lation deals with distribution functions and the momen- 
tum relaxation times are for ensemble averages only: The 
longer momentum relaxation time results in the minority 
electrons arriving earlier than the majority ones, which 
causes a characteristic shape of the spin dynamics with 
a negative spike in the spin density. When the major- 
ity electrons arrive, the signal of the minority carriers 
already falls off. Therefore the spin density n m changes 
sign and only afterwards shows a long positive tail. Im- 
portantly, the difference between the momentum relax- 
ation times can be changed from 1% to 100% leading to 
changes in the numerical values of the spin density, but 
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FIG. 5. Time evolution of the spin density of majority elec- 
trons n + (z = L, t) (a) and the total spin density n m (z = L, t) 
(b) for Au layer thicknesses of L = 25 nm (solid), 50 nm 
(dashed), and 100 nm (dotted). 



not in the qualitative behavior. Thus, as long as there 
is a difference between r+ and r_ , the resulting shape of 
the spin dynamics with the characteristic short negative 
dip and a long positive decay is extremely robust. 

We now discuss the results for the case in which the 
finite thickness L of the metal layer is modeled by a 
boundary condition for n s (z,t) at z = L. Figure [5^ a) 
shows the spin density n + (z = L,t) for majority elec- 
trons at the right boundary of the gold slab as a function 
of time for slab thicknesses L = 25, 50, and 100 nm, re- 
spectively.^ Several oscillations are visible for the 25 nm 
slab. For the 50 nm slab, only a second peak occurs be- 
fore the signal decays. This happens because, unlike the 
case in Fig. |4|a), the electrons are reflected at the right 
metal surface and continue to propagate/diffuse through 
the slab. Depending on the thickness of the layer, they 
may be reflected also at the left boundary, i.e., the Fe/Au 
interface in the experiment, and propagate again to the 
right metal boundary. This behavior continues as long as 
there is a ballistic component to the dynamics and van- 
ishes when the dynamics become diffusive. Eventually, 
the signal decays on the time scale of the spin relaxation 
time. In the 100 nm thick layer, in particular, the bal- 
listic character is further reduced, so that no additional 
peaks develop at all. The additional peaks, if present in 
the density dynamics n + of one spin channel, are also 
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z (nm) 

FIG. 6. Snapshots of the spin-current density J+{z,t) (a) 
and of the spin density n + (z, t) (b) as a function of z at times 
t = Qfs, 50 fs, 100 fs, and 150 fs for a 100 nm gold slab. 



visible in the total spin-density dynamics, n m , as shown 
in Fig.[5]jb). Depending on the thickness of the layer, the 
spin-signal can change sign more than once for the same 
reason as without reflection because the propagation and 
reflection dynamics occurs for both spins with a slightly 
different relaxation time, so that the maxima and min- 
ima of ri- are slightly shifted with respect to those of 
n+. 

Figure[6]shows the space dependence of the current and 
density distributions for "+" electrons for a slab of thick- 
ness L = 100 nm at different times. The snapshots of the 
current density, Fig. |6|a) show how the initial current 
pulse propagates through the gold slab, is attenuated, 
and changes sign after it reaches the right slab boundary. 
The sign change indicates that the propagation direction 
of the carriers is reversed, i.e., carriers moving ballisti- 
cally are reflected when they reach the boundary. This 
is another indicator that even in the case of the 100 nm 
gold slab, the ballistic character of the transport is still 
recognizable at the surface. The spin density n + (z,t), 
shown in Fig. |6^b), exhibits a wave front at first and a 
subsequent "smoothing" out due to the diffusive charac- 
ter of the dynamics, which becomes dominant after a few 
hundred femtoseconds. 

Finally, Fig. [7] shows the time evolution of the charge 
density n(z — L,t) at the right boundary of the gold 
slab for slab thicknesses L = 25 nm, 50 nm, 100 nm and 
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FIG. 7. Time evolution of the charge density n(z, t) for layer 
thicknesses of L = 25 nm (solid line), 50 nm (dashed) 100 nm 
(dotted), and 150 nm (dash-dotted). 



150 nm, respectively. Dependent on the thickness there 
are additional peaks due to multiple reflections on the 
boundaries and charge accumulation in the gold layer. 
The charge accumulates and reaches a steady state be- 
cause we neglect a possible flowing back of the charge- 
current into the iron layer. In an infinitely thick gold 
layer, the Gaussian current pulse, which mimics the ex- 
citation, propagates through the nonmagnetic gold while 
its maximum decays due to the spin-diffusion length. 
Thus, the space-resolved results for the spin and charge 
dynamics show that there is no spin and charge accumu- 
lation in an infinitely thick gold layer. 



C. Relation to experiment 

We have already described the idea of the experiment 
by Melnikov et al., and our calculation applies to the ba- 
sic setup of that experiment, as shown in Fig. [3| It is 
important to note that we do not calculate the electric 
fields that are generated by second-harmonic effects at 
the back surface and whose intensities are measured in 
the experiment. However, we follow the interpretation of 
the experimental results by Melnikov et al., and compare 
the "charge" signal with the charge density n and the 
"magnetic" signal with n m . An important experimen- 
tal finding is a fast rise and a subsequent slow decay of 
the charge signal accompanied by a magnetic signal that 
shows a characteristic sign change. Comparing the spin 
density n m at the right boundary for slab thicknesses 
50 nm and 100 nm, see Fig.|4^b) with the magnetic SHG 
signal shown in Figs. 2(c) and (d) of Ref. ITSIone finds a 
qualitative agreement. In particular, the time at which 
the sign change of the signal occurs and the decay at 
longer times are quite well described by the calculated 
result. Melnikov et al. observe that their magnetic sig- 
nal goes to zero at about 1.2 ps and therefore give roughly 
1 ps as an estimate of the spin-relaxation time (denoted 
as t^ u there) pertaining to their experiment. With the 
help of our model, we can improve this estimate consid- 
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erably, because our spin relaxation time Ti = 515 fs is a 
parameter of the calculation and its numerical value has 
been determined from transport measurements in the dif- 
fusive regime. The calculation therefore provides a link 
between the ultrafast, optically excited/detected dynam- 
ics, and "conventional" transport studies. Regarding the 
degree of agreement, we would like to add two remarks. 
First, the relative height of the extrema in the calculated 
spin-density dynamics is somewhat off, but it is possible 
to improve the agreement with regard to the peak height 
by making the excitation spin dependent, i.e., by having 
J_l_ 7^ J_ . This is likely the case in the experiment but 
including a spin dependent excitation would introduce an 
adjustable parameter, and we are here not interested in 
fitting the data, but in analyzing the physics contained 
in the macroscopic equations and Second, the 
calculation seems to lead to sharper extrema and smaller 
signals at longer times. Very likely, this is due to our ne- 
glect of energy relaxation processes because the energy 
dependence of the distributions underlying the macro- 
scopic equations is assumed to be essentially fixed and 
concentrated around the Fermi surface.^ The last point 
is also important when comparing the charge dynamics 
in our Fig. |4][a) and Figs. 2(a) and (b) of Melnikov et 
al. The area under the calculated curves in Fig. Qa) 
decreases strongly for larger values of L, whereas in the 
measurement it increases. In the calculation, this is a 
consequence of the spatial redistribution of the injected 
carrier density, as can be seen in Fig. [6jb) . In the exper- 
iment it is probably a result of the energy dependence of 
the nonlinear process that gives rise to the measured sig- 
nal, and an energy-redistribution of carriers during the 
dynamics. It is hard to say which is more important 
because a microscopic description of the SHG is quite in- 
volved.^ To have a better comparison of the qualitative 
behavior of charge dynamics for the different slab thick- 
nesses, we normalize the computed charge dynamics in 
Fig. [5J The normalization not only removes the differ- 
ence in the peak heights, but also leads to a behavior at 
long times for both slab thicknesses that is very similar 
to the experimental result. 

The qualitative agreement of the of the spin and 
charge dynamics with the signatures seen in the exper- 
iment reveals, in our opinion, a view of the experiment 
that is complementary to that given by Melnikov et al. 
There, the energy dependence of the excited carriers was 
stressed, along with the fact that one needs both elec- 
trons and holes with lifetimes that have to match up 
rather well to be able to separate carriers into those that 
behave mainly diffusively and those that behave mainly 
ballistically. Our simple macroscopic model that con- 
tains both a ballistic and diffusive transport properties 
and using well-established transport parameters suggests 
that the overall behavior seen in the experiment does not 
depend critically on details of the electronic lifetimes. 
Instead, qualitative signatures of the observed dynamics 
are likely quite generic for transport in thin metal slabs 
on ultrashort time scales, regardless of the microscopic 
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FIG. 8. Scaled dynamics of the density n+(z,t) of major- 
ity electrons at points z = 50 nm (dashed), 100 nm (dotted), 
150 nm (dash-dotted), and 300 nm (solid) of an infinitely thick 
layer. 



details of the excitation and relaxation dynamics. 

Finally, we compare the general trends seen in Fig. [8] 
with a pioneering measurement!^ of heat transport in 
thin gold films excited by a femtosecond laser pulse. This 
stretches the limits of the applicability of our calcula- 
tion because in that experiment intraband excitation was 
likely dominant, leading only to a heating of the conduc- 
tion electrons and not to an injection of carriers. How- 
ever, the macroscopic equations for heat transport have 
th e sam e mathematical properties as the ones used by 
us p5 | i6 | g0 ifoafc one can try to find indications of a qualita- 
tive agreement. In this case, we have to compare the nor- 
malized dynamics as shown in Fig. [8] because their Fig. 2 
also shows normalized reflectivity changes. Indeed, the 
overall shape of the reflectivity changes observed by Bror- 
son et al. match up well with our results. In particular, 
in Ref . [T3] it is found "that the rise time of the signal in- 
creases slightly with increasing thickness." This trend is 
clearly observed in our calculated results and is explained 
by a change of dominantly ballistic dynamics in thinner 
samples to dominantly diffusive dynamics in thicker sam- 
ples, here 300 nm. Importantly, this transition comes out 
of the calculation by a combination of ballistic and diffu- 
sive aspects: Even in thicker samples where the diffusive 
dynamics dominates the rise time of the signal, there is 
still an important component of the result that can only 
be explained by ballistic dynamics, namely the delayed 
onset of the reflectivity response. Thus it is hard to ex- 
plain the observed rise times in terms of either ballistic or 
diffusive transport, and Brorson et al. noticed this incon- 
sistency in the interpretation of their data. An analysis in 
terms of wave-diffusion dynamics removes the apparent 
contradiction between diffusive and ballistic properties of 
transport. 



V. CONCLUSIONS 



This paper is devoted to an analysis of the wave- 
diffusion equations, Eqs. ([!]) and (|2|), which provide a 
macrosopic description of carrier transport due to optical 
or electrical injection. We showed that the wave charac- 
ter in Eqs. ([lj and ([2| is related to ballistic transport by 
analyzing the behavior of an analytical solution to these 
equations in one dimension. Scattering processes change 
the behavior from ballistic at short times to diffusive on 
a time scale larger than the momentum relaxation time. 
Importantly, we find that our simple, macroscopic equa- 
tions work not only in the purely ballistic and diffusive 
limits, but that they can also applied on time and length 
scales, where spin and charge transport are neither com- 
pletely ballistic nor diffusive. Using the wave-diffusion 
equations, we model the conditions realized in a recent 
experiment, in which spin-polarized electrons are excited 
in a gold film by optical pumping of an adjacent ferro- 
magnetic layer and probed optically on the other side 
of the gold film. If we assume different relaxation times 
for spin-up and spin-down electrons, our simple macro- 
scopic model can reproduce qualitatively important sig- 



natures seen in the experiment, which underscores the 
importance of the transition between ballistic and dif- 
fusive aspects of transport for an understanding of the 
observed dynamics. Our model uses established macro- 
scopic transport parameters and the dependence of our 
results on the one additional parameter that we need to 
introduce, namely the difference of the relaxation times, 
is very weak. Our results are thus robust, which indicates 
that the observed behavior is quite generic for ultrafast 
transport dynamics in thin metal films. Last, we com- 
pared the overall behavior of our calculated results to a 
reflectivity measurement on gold films, and suggest that a 
unified description of ballistic and diffusive transport can 
clarify the behavior of electronic transport at ultrashort 
timescales that may seem unintuitive when one tries to 
classify transport as either ballistic or diffusive at differ- 
ent times. 



ACKNOWLEDGMENTS 

We acknowledge support by the state of Rhcinland- 
Pfalz through the MATCOR program and the MAINZ 
graduate school, as well as by the Beijing Natural Science 
Foundation (no. 1112007). 



* hcsch@physik.uni-kl.de 

1 I. Zutic, J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 
76, 323 (2004). 

2 J. Fabian, A. Matos-Abiague, C. Ertler, P. Stano, and 
I. Zutic, Acta Phys. Slov. 57, 565 (2007). 

3 A. Brataas, G. E. W. Bauer, and P. J. Kelly, Phys. Rep. 
427, 157 (2006). 

4 M. Johnson and R. H. Silsbee, Phys. Rev. Lett. 55, 1790 
(1985). 

5 F. J. Jedema, A. T. Filip, and B. J. van Wees, Nature 410, 
345 (2001). 

6 T. Valet and A. Fert, Phys. Rev. B 48, 7099 (1993). 

7 S. Zhang and P. M. Levy, Phys. Rev. B 65, 052409 (2002). 

8 E. I. Rashba, Appl. Phys. Lett. 80, 2329 (2002). 

9 J. Zhang and P. M. Levy, Phys. Rev. B 71, 184417 (2005). 

L. Cywihski, H. Dery, and L. J. Sham, Appl. Phys. 
Lett. 89, 042105 (2006). 

1 Y.-H. Zhu, B. Hillebrands, and H. C. Schneider, Phys. Rev. 
B 78, 054429 (2008). 

2 Y.-H. Zhu, B. Hillebrands, and H. C. Schneider, Phys. Rev. 
B 79, 214412 (2009). 

3 A. Melnikov, I. Razdolski, T. O. Wehling, E. T. Papaioan- 
nou, V. Roddates, P. Fumagalli, O. Aktsipetrov, A. I. 



Lichtenstein, and U. Bovensiepen, Phys. Rev. Lett. 107, 
076601 (2011). 

14 S. D. Brorson, J. G. Fujimoto, and E. P. Ippen, Phys. Rev. 
Lett. 59, 1962 (1987). 

15 D. D. Joseph and L. Preziosi, Rev. Mod. Phys. 61, 41 
(1989). 

18 M. Chester, Phys. Rev. 131, 2013 (1963). 

17 F. Steininger, A. Knorr, P. Thomas, and S. W. Koch, Z. 
Phys. B 103, 45 (1997). 

18 G. H. Weiss, Physica A 311, 381 (2002). 

19 M. Battiato, K. Carva, and P. M. Oppeneer, Phys. Rev. 
Lett. 105, 027203 (2010). 

20 V. V. Kruglyak, R. J. Hicken, M. Ali, B. J. Hickey, A. 
T. G. Pym, and B. K. Tanner, Phys. Rev. B 71, 233104 
(2005). 

21 Wen-C. Chiang, C. Ritz, K. Eid, R. Loloee, W. P. Pratt, 
and J. Bass, Phys. Rev. B 69, 184405 (2004). 

22 The numerical values of n s (z,t) are twice as large as the 
ones in Fig. [4] because there is no carrier density "lost" to 
propagation outside the slab. 

23 G. Lefkidis and W. Hiibner, Phys. Rev. Lett. 95, 077401 
(2005). 



